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Summary. Early-type galaxies exhibit a wealth of photometric and dynamical 
structures. These signatures are fossil records of their formation and evolution pro- 
cesses. In order to examine these structures in detail, we build models aimed at 
reproducing the observed photometry and kinematics. The developed method is a 
generalization of the one introduced by [l^ ]. consisting in an N-body representation, 
in which the weights of the particles are changing with time. Our code is adapted for 
integral-field spectroscopic data, and is able to reproduce the photometric as well as 
stellar kinematic data of observed galaxies. We apply this technique on SAURON 
data of early-type galaxies, and present preliminary results on NGC 3377. 



1 Introduction 

Goal Elliptical galaxies are presumed to be the product of galaxy merg- 
ing, where different mass ratios result in different classes of elliptical galaxies 
(disky, boxy. . . ) M. This seems consistent with the hierarchical scenario of 
galaxy formation [9(, although it is clear that more standard collapse pro- 
cesses and secular evolution contribute to the shaping of galaxies. There are 
many different and often complementary approaches to constrain these sce- 
narios. The one we adopted is to study the structures of nearby elliptical 
galaxies in detail, in order to probe the signatures of their evolution and 
to be able to trace their formation history. This fossil search requires data 
both on the chemical and dynamical status of the stellar components, and if 
present of the gaseous system, both being provided by spectroscopy. 

3D spectroscopy Spectroscopy can indeed deliver information on the stars' 
dynamics and chemistry via the analysis of their absorption lines. Two- 
dimensional spatial coverage provided by integral-field spectroscopy is critical 
to homogeneously sample an extended target, while more traditional long- 
slit spectroscopy generally restricts our view to a few a priori fixed axes. A 
qualitative assessment of the obtained maps is clearly not sufficient to fully 
address the issues mentioned above. We need to go further by modelling the 
observed galaxies using constraints retrieved from 3D spectroscopy, such as 
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luminosity, stellar velocity distribution and stellar populations. As empha- 
sized by [lj, 3D spectroscopy seems to be essential to properly constrain the 
dynamics of the galaxy under scrutiny. 

Existing models One of the key goal of galaxy modelling is to retrieve the 
full distribution function (DF), i.e. the density of stars in phase-space. There 
are already several known methods to address this problem, each of them 
having both advantages and drawbacks. We can distinguish different classes 
of techniques : 

DF-based method, where models are often restricted to simple geometries. 
Moment-based method, which consists in solving Boltzmann and Poisson 
equations via the use of a closed system of relations (Jeans equations). 
The main issue here is that the final DF may not be p ositive everywhere. 
Orbit-based method, or Schwarzschild modelling [10(. Libraries of orbits 
are built within a fixed potential, and each orbit is weighted as to repro- 
duce the observed galaxy (photometry and kinematics). 
Particle-based method, where the particles represent groups of stars. One 
has to guess the right initial conditions which will evolve in a configuration 
resembling the chosen galaxy. 

2 Method 

The method we wish to present here is an hybrid scheme between the 
Schwarzschild and N-body techniques. It has been first proposed and de- 
veloped by [13] and consists in an N-body realization where the weight of 
each particle is gradually changed in order to fit the observables. 

Previous applications This method has been tested by [T3| on fake galaxy 
models, adjusting the photometry alone. It has been recently applied on the 
Milky Way by [2| , fitting the photometry of a frozen snapshot extracted from 
a full N-body simulatiorQ- 

Algorithm The algorithm we use is illustrated in Fig. Q] Particles start from 
(N-body) initial conditions and are integrated along their orbits. Observables 
are used as reference input to the modelling, particles being projected as to 
mimic these observables. A weight prescription is derived from the difference 
between the model and the data, this prescription being then applied on 
each particle accordingly. This procedure (integration, projection, compari- 
son, weight changing) is repeated until the weights converge. 



1 An extended independent implementation of this algorithm has been recently 
presented by de Lorenzi et al. 2007: see their paper for details 
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Fig. 1. Algorithm of the Syer & Tremaine method. 



Initial conditions The initial conditions have to properly sample the galaxy, 
ie in position/velocity space or using integrals of motions. This is not an easy 
task when we do not a priori know the internal structure of the galaxy we 
wish to model. If we want to model real galaxies, we have to find appropri- 
ate initial conditions, without a priori knowing its DF. We thus first build 
an initial mass distribution, by using the known photometry and assuming 
a constant mass-to-light ratio. This is achieved here via the Multi Gaussian 
Expansion (MGE) proposed by [7j, which allows us to decompose the mass 
distribution into a sum of tri-dimensional Gaussians, and to derive the cor- 
responding gravitational potential analytically (f|. Jeans equations are then 
used to obtain a first guess of the dynamics and set up initial conditions. 



Integration In a first version of the code, we keep the potential steady. The 
code has been developed in a very modular and flexible fashion, such as to 
allow the easy implementation of a self-consistency module : this will be de- 
scribed in a forthcoming paper. As we need to evaluate the observables of 
all particles at the same time, we chose a second-order synchronized leapfrog 
scheme, where positions and velocities can be evaluated simultaneously. The 
leapfrog technique is perhaps not the best choice in terms of accuracy, but was 
found sufficient for this first development step. Other schemes (e.g. Runge- 
Kutta) can be easily implemented. Adaptive time-steps have been included, 
in order to optimally sample all the orbits. An independent integration of 
each particle would result in a very inefficient algorithm. We thus chose to 



follow the algorithm described in the N-body code GADGET 11], in which 
particles advance in bunches so that they are always distributed in a tight 
time range around the current time. 



Observables Our code has been designed around two new main items. First, 
we do not restrict ourselves to fit the photometry : the code also allows the 
fitting of the kinematics, via the use of Line of Sight Velocity Distributions 
(LOSVDs). Second, the code has been developed as to permit the use of 
3D spectroscopic data, including spatial adaptive binning (Voronoi bins, see 
0]). We have thus used the data obtained with the SAURON spectrograph 
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(WHT, Canary Islands) for its survey of early- type galaxies 0,0. 



Prescription The heart of the code is the prescription developed by [12| . 
which rules the weight changing of each particle so that the cumulated ob- 
servables reproduce the observations : 



where i describes the particles, j the observables, w the weights, the 
contribution of the particle i to the observable j, and Aj = yj(t)/Yj — 1 the 
contribution of all the particles to this same observable, e is the strength of 
the weight changing, Wi avoids the weights to be negative when Wi approaches 
0. The input parameters are the following : 

- N is the number of particles, 

dtime is the time during which the observables are averaged, 
e is the strength parameter mentioned in Eq. [IJ 

a is the smoothing parameter of the observables, used in the equation : 



so that Aj (t) replaces Aj (t) in Eq. [TJ 

The choice of these parameters is crucial, each galaxy requiring an appropri- 
ate set. 

3 Results 

A preliminary model has been obtained on the nearby disky galaxy NGC 3377. 
In Fig. [2] we illustrate the fit of the photometry and LOSVDs obtained with 
the SAURON instrument on this galaxy, for a given mass (MGE) distribution. 
The maps are globally well reproduced (I, a, h^). Some remaining discrep- 
ancy in the velocity field is probably due to an observed asymmetry in the 
kinematics of this galaxy (tilt between the photometric and kinematic minor- 
axis) . The relatively high /14 observed level is not reached by the model : this 
may be partly related to a template mismatch effect which mostly affects the 
even moments of the LOSVD. 

4 Conclusion and perspectives 

Although the results shown here are very preliminary, they illustrate the 
possibility of fitting the kinematics, which is an important step in the devel- 
opment of such a code. It has also been designed to adjust complex data such 
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as the one provided by 3D spectrographs. A few improvements have been or 
are now being implemented : 

Self-consistency : it is technically easy to replace the current integration 
module with a self-consistent scheme. The full implications of a continuous 
weight changing remains however unclear and have to be examined. 
Addition of a central black hole : this would allow the use of higher spatial 
resolution data of nearby galaxies. 

Stellar populations : this new class of observables would allow us to link 
the dynamics and chemistry of the stellar component. 
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